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Abstract 

Closed loop solitons in a plane, whose curvatures obey the modified Korteweg-de Vries equation, 
^5 were investigated. It was shown that their tangential vectors are expressed by ratio of Weierstrass 
^ gigma functions for genus one case and ratio of Baker's sigma functions for the genus two case. 
T^lj- This study is closely related to classical and quantized elastica problems. 
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1. Introduction 



This article is on loop solitons of genera one and two. We give them as ratios of sigma functions; 
for the case of genus one, the sigma functions are Weierstrass sigma functions while for the case of 
genus two, they are Baker's hyperelliptic sigma functions. 

First wc mention background of the study, which is related to elastica problem. Elastica problem 
was proposed by James Bernoulli in ending of seventeenth century and Euler and his nephew Daniel 
Bernoulli investigated the problem [T, L, W]. The elastica problem is to determine all allowed shapes 
of a thin elastic beam in a plane. The allowed shape is realized as a curve with the local minimum 
of the energy 

E = Jdsk^, (1-1) 

where s is its arclength and k is its curvature. For its tangential angle (p and g := ex.p{^/—l(p), the 
integrand is given by 

k'^ds = -ig-'dg){*g-'dg) = id,<f>fds, (1-2) 

where * is the Hodge star, a map from one-form to zero-form, and dg := d/ds. Thus the elastica 
problem is the oldest harmonic map (sigma-model) problem from a circle or a real line M to 
U(l). This problem was essentially solved by Euler using numerical computations [T,L]. 

For a given elastica in two dimensional euclidean space, it is characterized by the affine coordinate 

(Xi(s),X2(s)) or Z{s) = X'^{s)+^/^X'^{s). Then its curvature k is given by k := "^/^^s logd^Z. 
The shape locally minimizing the energy (1-1) is given by the differential equation [Mai, Ma2], 

Cd,k + ^k'^d.k + dlk = 0. (1-3) 

This is known as the static modified Korteweg-de Vries (SMKdV) equation. 

As I showed in ref . [Ma2] , when we deal with a non-stretching closed elastica in thermal bath or 
investigate the partition function of temperature 1//3, 



Z = VZexp 



{-13 J dsk^), (1-4) 



we also need information of curves which have the excited energy. We will call such curves "quan- 
tized elastica" following [Ma4]. This problem is associated with large polymer problem like DNA 
in biology and polymer physics. As shown in [Ma2,Ma3,Ma4], the quantized elastica is classified 
by the modified Korteweg-de Vries (MKdV) equation, 

Cdtk + ^k'^dsk + 5f A; = 0, (1-5) 

where ds '■= d/ds, dt := d/dt and t is a deformation parameter. The orbits obeying the MKdV 
equation preserve local length and their first integral (1-1), i.e., [dgjOt] = and dtE = 0. Since 
the MKdV equation is an initial value problem, for a given smooth curve, there exists, at least, 
an MKdV orbit including the curve; wc can classify quantized elastica using the orbits. Hence the 
partition function (1-4) is reduced to an integration. 



J dn{E)exp{-(3E), (1-6) 



where dQ{E) is the density of states. The evaluation of dO,{E) means the determination of the 
orbits of the MKdV equation whose first integral is E. In other words, in order to evaluate (1- 
5), we must explicitly construct the periodic solutions of MKdV equation, which are expressed 
by hyperelliptic functions. The hyperelliptic function is characterized by genus. As mentioned in 
[Ma2], orbits with small genus is expected to largely contributes to the partition functions (1-4). 
Thus as an attempt, we will investigate the classical and quantized elastica of genera one and two 
respectively. 

Since the curve whose curvature obeys the MKdV equation was called loop soliton [WIK, I], 
our investigation is, just, a study of closed loop solitons. Here we will also note that the time 
development of classical elastica does not obey the MKdV equation except genus one solution case 
[Mai]. (From physical point of view, investigation of a loop soliton with higher genus is somewhat 
nonsense in pure kinematic study.) 

In this article, we will show that the tangential vectors of closed loop solitons with genera one and 
two are expressed by ratio of Weierstrass sigma functions in §3 and ratio of Baker's sigma functions 
§4 respectively. This result is very interesting as we will discuss in §5. The sigma function, which 
is well-tuned theta function, plays very important roles in the studies of elliptic and hyperelliptic 
functions [Bl, B2, B3, BEL, G, Ma5, 01, 02, 03]. 

Further our result has an effect on historical interpretation of history of science. According to 
Truesdell and Love [T,L], Euler did illustrate the shape of loop soliton with genus one, our result 
means that Euler essentially found the elliptic functions and sigma functions even though he dealt 
only with elliptic integrations. In fact, the tangential vector dsZ{s) is periodic function of s. I 
believe that Euler noticed its properties and thus after he read Fagnano's article, his mind felt the 
addition formula [W]. 

The construction of closed loop soliton with genus two is based upon Baker's hyperelliptic 
function theory [B1,B2,B3]. Recently Baker's approach is reevaluated from various viewpoints 
[BEL, G, Ma5, 01, 02, 03]. Buchstaber, Enolskii and Leykin has been actively studying the 
Baker's theory from soliton theory and developing it. Grant and Onishi found interesting relations 
in hyperelliptic function in terms of Baker's theory from point of view of number theory. As on 
Baker's theory of the hyperelliptic function, we start from a hyperelliptic curve and find an explicit 
function form without ambiguous parameter, it is very effective. 



Let us consider a smooth immersion of a circle S into the two dimensional euclidean space 
C. The immersed curve C is characterized by the affine coordinate {X^{s),X'^{s)). Here s is 
a parameter of and is, now, chosen as the arclength so that ds^ = (dX^)^ + {dX^y. We will 
also use the complex expression, 



2. Basic Properties of Loop Soliton 



Z{s) :=X\s) + V^X\ 



(2-1) 



Then |5sZ(s)| = 1 and the curvature of C is given by 



kis) 



1 



dslogdsZ{s). 



(2-2) 



3 



Definition 2-1. 

A one parameter family of curves {Ct} for real parameter t G M «s called a loop soliton, if its 
curvature obeys the MKdV equation; for q := k/2, 

dtq + 6q^d,q + = 0. (2-3) 

Here we will describe the Miura map for later convenience. 

Proposition 2-2 [DJ]. 

For the solutions p± of the KdV equation, 

dtp± + 6p±dsp± + dlp± = 0, (2-4) 

if we find a quantity q satisfying p± = q^ ziz dgq over Ct and for all t, q obeys the MKdV equation 
(2-3). 



3. Loop Soliton with Genus One and Weierstrass Sigma Function 



In this section, we will deal with a loop soliton with genus one. In other words, we are consider 
the case oit = s/C, i.e., the SMKdV equation, 

Cd,q + &q^d,q + g = 0. 



First we will set up the Weierstrass p-functions and cr-f unctions. 

Definition 3-1 [WW]. 

For an elliptic curve, 

= A{x - ei)(x - e2){x - 63), 

(1) We will define integrals u, Ua, rja (a = 1,2,3^, 

[^""'^^ dx /•("-°) dx /•(^-°) xdx 
u:= — ^ oja= — , ria= 

^(0,0) y J(o,o) y J{o,o) y 

(2) The elliptic theta function of t := ujz/loi is defined by 
9^(z) := ^exp (iV^tt 
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(3) The Weierstrass a function is defined by 

(t{u) := exp(?7i'u^/(2a;i)) 

(4) We also introduce the other a-functions, 



(V(2^i)) 



du0i{u/(2ui))\u=o' 



aa(w) := exp(-77„w) — — , a = 1,2, 3. 

o"('^aj 

(5) The Weierstrass p-function is defined by 

p := -dlloga{u). 



(3-1) 



(3-2) 
(3-3) 

(3-4) 

(3-5) 

(3-6) 
(3-7) 
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Theorem 3-2. 

The shape of loop soliton with genus one is given by 



^ ^ J \ais- W3/2 + S) J ' ^ ^ 

where S is a constant parameter. 



Later part in this section devotes the proof of this theorem but first we wih summary the 
properties of Weierstrass (elhptic) p and a functions. Let ""' be derivative in u. 

Proposition 3-3 [WW]. 

(1) 

p'iuf = 4p^ + g3p + gi 

= 4(p-ei)(p-e2)(p-e3), (3-9) 

where p{(jJa) = o,n-d ei + 62 + 63 = 0. 

(2) 

12p{u)dMu) + dlpiu) = 0. (3-10) 
Further we will note the addition relations. 



Proposition 3-4 [WW]. 

(1) 



(2) 



(3) 



(4) 



o'iz + u)a(z — u) 



(u + .j-pW^-k lM-f/"' . (3-12) 
2 p{z) - p[u) 



Here (3-12) and (3-14) is directly derived from (3-11). (3-13) is slightly difficult but using 
proposition 3-3, it can be proved. 



Prom the definition of loop soliton, the theorem 3-2 is reduced to the next lemma. 
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Lemma 3-5. 

For p := p{u) — 63, ^i{u) := ^ du logp obeys the MKdV equation, 

duU + Gfi^duH + din = 0. (3-15) 

Proof. Prom proposition 3-4, 

p{u + U3) - p{u) = --j==duH, (3-16) 

p{u + u;3) + p{u)+e3 = -iJ,'^. (3-17) 
Thus for V := u — CO3/2, we have a relation, 

-2p(t; ± UJ3/2) - 63 = /x^(t;) ± -\/^5^/x(t;). 

Prom (3-10), p = — 2p(u) — 63 satisfies the equation, 

6e3duP + GpduP + dip = 0. (3-18) 

Using the proposition 2-2, this lemma is proved. I 

Hence the theorem 3-2 is also satisfied. 
Remark 3-6. 

We note that (3-8) is a necessary condition. If one tunes S in (3-8) and the g2 and 53 in (3-2) 
to find a real line s in the complex plane {u} = C so that it satisfies the reality condition \dsZ\ = 1 
and the closed condition Z{s) = Z(s + 2loi), (3-8) becomes the closed loop soliton. 



i4. Loop Soliton with Genus Two and Baker's Sigma Function 



In this section, we will deal with a hyperelliptic solution of the loop soliton whose genus is two. 
We will redefine the quantities cr, p, p, uj, r], t in the previous section. These quantities are altered 
to genus two versions of corresponding quantities in the previous section. 

Definition 4-1 [Bl, B2, BEL Chapter 2, 61 p.384-386, 02]. 

For a hyperelliptic curve of genus two, 

= fix) 

= Ao + AlX + A2X^ + ••• + A5X^ (4-1) 

where A5 = 1 and Xj 's are complex numbers, we will give definitions as follows. 
(1) We choose ar, Cr and c {r = 1,2) for 

f{x) = P{x)Q{x), 



P{x) = {x — ai){x — 02), Q{x) = {x - ci){x — C2){x — c), 
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(4-2) 



so that curves (a^, c^) or (c, oo) does not intersect each other and Re{ar) < Re{cr), Re{ai) < 
Re{a2), Re{c2) < Re{c). 
(2) Let us denote the homology of the hyperelliptic curve by 



Hi(Xg,Z) = Zai e Z/3i e Za2 Z/^a, 

where these intersections are given as [ai,aj] = 0, [Pi,(3j] = and [ai,(3j] = Sij. 

(3) The unnormalized differentials of first kind are defined by, 

dx xdx 
dui := — , du2 := — — . 
2y 2y 

(4) The unnormalized differentials of second kind are defined by, 

dri := — (A3X + 2X4X^ + SX^x^) dx, dr2 '■= — X^x^dx. 
2y ^ 2y 

(5) The unnormalized period matrices are defined by, 

dui 



2w' :-- 



/ dui 


/ dui 






.7 a 2 




/ dU2 


/ dU2 




J ai 


J 012 





/ du2 I 

JP-L J 132 



du^ 



duo 



u> :- 



(6) The normalized period matrices are given by, 



/-I // 
T ■.= 10 (J , U) :-- 



(7) The complete hyperelliptic integral of the second kinds is given as 

dr 



2'n' :-- 



dri I o'/'i 

CKi J 0:2 



dr2 



2 / dr2 

Oil CK2 



dri / d^i 

111 J 132 

dr2 / dr2 

01 Jl32 



(8) We will defined the Riemann theta function over characterized by A, 



{z) :- 



{z-T) 



exp 



27rV-l {i^\n + a)T{n + a) + \n + a){z + b) 



for 2- dimensional vectors a and b. 
(9) We will introduce constant parameters, [B2 p. 343] 



Xr 



(V-l/'(a.)/4)V2^ 



r = 1,2. 



(4-3) 



(4-4) 



(4-5) 



(4-6) 



(4-7) 



(4-8) 



(4-9) 



(4-10) 



We will note that these contours in the integral are, for example, given in p. 3. 83 in [M]. Thus 
above values can be, in principle, computed in terms of numerical method for a given y^ = f{x). 
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Proposition 4-2 [Bl, B2 p.316-p.335, BEL p.25, M p.3.80-84, 02]. 

The Riemannian constant if G C^, G are given by 



/.(ai,0) l'l,a2,U) 

K = u'-\ du+ du) = 5' + 6"T, 



(a2,0) 



(4-11) 



where 



U>2 



Aai,0) 

:= / du = u'S[+u"S'^, 

•J oo 
/•(a2,0) 

:= / du = (j'S'2+cj"S'^, 

J oo 



6[:- 



- 1 - 




r ^ 1 




- 1 - 


2 




2 


, 5' = 






1 


1 


.0. 




- 2- 




-2- 



- 1 - 








r ^ 1 


2 








2 






1 




1 


. 0. 




-2- 




-2- 



(4-12) 



(4-13) 



Definition 4-3 [Bl p.286, B2 p.336, p.353, p.370, BEL p.32, p.35, Ol p.386-7, 02]. 

We will define the coordinate (1x1,^2) in for points (xi, yi) and (x2, 1/2) of the curve = f{x), 

Uj := / duj + / duj. (4-14) 



(1) Using the coordinate uj, sigma functions, which are homomorphic functions over 'C? , are 
defined by, 



a{u) = a{u;(j) := 7exp(— ^ * wq'u' ^uj-d 



6" 
S' 



ii2wr'u;T), 



(4-15) 



ar{u) = Ur{u\u}) := Ar-exp(— - * urji'u}' ^u)'& 



5'; 



((2a;')-'u;T), 



for a = 0,1, 2, 



7 ^0 = 1; 



where j is a fixed constant, Sq = Sq := 
(2) In terms of a function, p-functions over the hyperelliptic curve are given by 



\oga{u), 
log ar{u). 



(4-16) 



(4-17) 
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Theorem 4-4. 

The shape of loop soliton with genus two is given by 



J \ao{u{s,t)+Ui+u;2/2 + 5)J ^ ' 



1 /-4t 
/or 02 = —0-^4? where S is a constant parameter and u{s,t) := I 

o V ^ 



Before proving this theorem, we will also review the properties of Baker's p and cr-functions. 

Proposition 4-5 [B3 p.155-6, BEL p.56-58, G p.103-4, p.ll6, 6l p.388,62. 

Let us express p^up dpnv{u)/dup and p^vpx := d^pfj,^{u)/dundui,. Then hyperelliptic p- 
functions obey the relations 

{H - 1) P2222 - 6p22 = 2A3A5 + 4A4P22 + 4A5P2I, 

{H - 2) P2221 - 6p22p2i = 4A4P21 - 2A5P11, 

{H - 3) P2211 - 4p2i - 2p22pii = 2A3P21, 

{H - 4) P2111 - 6p2ipii = -4A0A5 - 2A1P22 + 4A2P21, 

{H - 5) pun - Gpji = -8A0A4 + 2A1A3 - I2A0P22 + 4Aip2i + 4A2pii. 



(I - 1) P222 = 4(p22 + P12P22 + A4P22 + Pll + -^3P22 + A2), 

(/ - 2) P222P22I = 4(pi2p22 - ■^iPllp22 - P?2 

+ A3P12 - Ai) + A3P12 + A4P12P22), 
(/ - 3) P221 = 4 (P11P22 - (PllP22 - P?2 + ^3Pl2 - Ai)p22 " Pllpl2 

+ A4P11P22 + A3P12P22 - A4(piip22 - P?2 + A3P12 - Ai) 

+A4A3P12 - A1P22 - A1A4 + Ao) . 

We also have the additional formulae for them 
Proposition 4-6 [B2 p.372, p.381, BEL p.ll2]. 

(1) 

a(u + v)a(u — v) / ^ , ^ / x / x / x / \ /. -,^\ 

. ^o . ^o = P22 (u) P21 (v) - P21 (u) P22 [v) - pii {u) + pii [v) . (4-19) 

(2) 

= T^M^i-) + P^2^(-)«2 - al). (4-20) 
ao{u) J Q{ar) 

We note that (4-19) corresponds to (3-11) and (4-20) is a genus two version of (3-14). The 
similar formulae of (3-13) and (3-14) are studied in [G p. 114-6] but they are not so effective for this 
purpose. For special value (theta characteristics), we give genus two versions (4-24) and (4-26) of 
(3-13) and (3-14) in proof of lemma 3-5. 
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Proposition 4-7. 



Proof. Using 9{z) := 



"0 



(z), the (7-functions are explicitly written by 



X t?((2a;'"^M + 5'+T5"), 



(4-21) 



ar{u) = Xr exp ( —-^ur]'u}' + 27r\/^ 



^*W+*<5;'((2a;')-^ + <5;)]) 



(4-22) 



Thus 



log(f7^(n)/(7(n)) - log(i9((2a;' \ + 6^ + T5';)/^ii2uj' \ + 6'+t5")) 



is a first order polynomial of u. Thus the derivative in the definitions of p leave only the difference 
of phases. Noting proposition 4-2 and K = (2a;')~^(a;i + UJ2), we obtain the identities, 

u = u — uji — UJ2 + 'iw'K, 
u + 2a;' ((5; + t5'1) = u-U2 + 2w'K, 



u + 2a;' ((^2 + T62) = u-ui + 2w'K. 

Hence (4-21) is proved. I 

Prom the definition 2-1, the theorem 4-4 is reduced to next lemma. 
Lemma 4-8. 

For a2 = —-X4, aiu) := , 9„„ log \ -^-tX\ obeys the MKdV equation, 

3 \<^o(w)/ 

Proof. By noting (4-20) and introducing 

p{u) := p^i^(^x) + p^2^(u)a2 - aj, 

fj,{s) is given by 



1 du2p{u) 



(4-23) 



2^/^ p{u) 

By taking the logarithm and derivative in M2, the addition formula (4-20) is reduced to 



(4-24) 

(4-25) 
(4-26) 



p22(« + (^2) - P22{u) = ^/-idu^niu + LUI+ UJ2). 
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(4-27) 



Thus we must prove the relation 

p22(li + ^2) + P22{u) + A4 + 02 = -tJL{u + Ui + LO-^)^ . (4-28) 

If (4-28) accomplishes, (4-27) and (4-28) give the relations that for v ■.= u — uj^l'^i 

-2p22{v ± ^2/'^) - A4 - 02 = IJ?{V + Wl + OO2) ± du^Kv + UJI+ UJ2). (4-29) 

Using (H-1) in the proposition 4-5, we can show that q := —2p22{v ± (^2/2) — A4 — 02 obeys 

-4du, q + 6qd^, q + d^^q = 0, (4-30) 

if 02 = — A4/3. From the proposition 2-2, if (4-28) is true, lemma 4-8 will be proved. 

Thus we will concentrate (4-28). We add (4-27) to 2p(u) -|- A4 -1-02 and then we have the relation 

p22{u + (^2) + P22('") + A4 + 02 = -n{u + UJi -|- ^2)^ + , (4-31) 

4p(« + UJI+ UJ2) 

where, 

A(u) := %p22{u)p(uf - 2{dl^p{u))p{u) + (4A4 + 03 + l){du,P2{u)f ■ (4-32) 

We use the relations in proposition 4-5 and then direct computation shows that A vanishes even 
though it is very tedious. Accordingly the lemma is proved I 



The theorem 4-4 is also proved and the remark 3-6 is also effective in this case. 



§5. Discussion 



We showed that the loop soliton is directly connected with a-functions. In ref.[Mal], 1 partially 
showed the fact but 1 could not reach the concrete function form and the difference between 
numerator and denominator in the ratio of a or (r) function. This study determines that the 
difference comes from the half of period or so-called theta characteristics [B2, M]. Even though we 
could not find a solution with higher genus, it is not difficult to conjecture the function form of 
loop soliton in terms of the theta characteristics. 

As we used the definition of loop soliton as a relation of MKdV equation, our study is just that 
of MKdV equation. In other words, we have concretely constructed the function forms of the genus 
one and the genus two periodic solutions of MKdV equation in the section 3 and 4. In ordinary 
study of elliptic function solution of MKdV equation [Mai], one uses Jacobi elliptic functions, 
sn, cn, and dn functions. Our solution (3-15) becomes the Jacobi elliptic function by means of 
Landen transformation, Jacobi transformation, Gauss transformation and so on. The relation 
between Weierstrass and Jacobi elliptic functions is derived from (3-14) [WW]. p{u) — plays 
very important roles; 



(\ / ei-eg p{u)-ei p{u) - 62 , . 

sn{z) = ^—— , cn{z) = \—— , dn{z) = \ —— , 5-1) 
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where z = (61 — 63)^/^1*. Since sn function is connected with 1/sn by phase shift and associated with 
dn and cn functions by algebraic equations, the Jaobi elUptic function could roughly be considered 
as p{u) — 63. Similarly we might regard ^/p in (4-25) as a Jacobi-type hyperelliptic function due 
to (4-20); ^fpjoa is essentially p(ii) — 63 when we set du^a = and u := U2- 

We should also notice the square in (3-14) and (4-20). As described in [Mai and reference 
therein], it is closely related to fermion. In fact, := (72 /cr ~ \/dsZ in (3-14) ( -0 := (72 /co iii 
(4-20) ) obeys the Dirac equation, 

(/I -aj (v^v) 

The square root of dgZ means double covering of the tangent bundle. (5-2) is related to Schwarz 
derivative and moduli space as mentioned in [Mai]. We recognize that there is still open problem 
what is the orign of the square root. 

Finally we will comment upon the closed condition mentioned in the remark 3-6. We could not 
touch the condition. In general, an integration of periodic function is not periodic. However if 
the function satisfies some conditions, its integral is also periodic. The closed condition is a such 
condition. It suppresses two degree of freedom in the parameter spaces. For the study of quantized 
elastica problem, it is very important to determine the relations between the conditions and the 
coefficients of the algebraic curve \j in (4-1). In future, we wish to determine them. 
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